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ABSTRACT 

We report on the Chandra X-ray Observatory and Hubble Space Telescope 
observation of two low core density globular clusters, NGC6144 and E3. By 
comparing the number of X-ray sources inside the half-mass radius to those 
outside, we found 6 X-ray sources within the half-mass radius of NGC6144, among 
which 4 are expected to be background sources; 3 X-ray sources are also found 
within the half-mass radius of E3, of which 3 is expected to be background 
source. Therefore, we cannot exclude that all our sources are background sources. 
However, combining the results from X-ray and optical observations, we found 
that 1-2 sources in NGC6144 and 1 source in E3 are likely to be cataclysmic 
variables and that 1 source in NGC6144 is an active binary, based on the X- 
ray and optical properties. The number of faint X-ray sources in NGC6144 and 
E3 found with Chandra and HST is higher than a prediction based on collision 
frequency, but is closer to that based on mass. Our observations strongly suggest 
that the compact binary systems in NGC6144 and E3 are primordial in origin. 
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novae, cataclysmic variables — X-rays: binaries 

1. Introduction 

The number density of X-ray sources in globular clusters with luminosities greater than 
L x = 10 36 erg s _1 is 100 to 1000 times higher than that of the Galactic disk (Giacconi 
et al. 1972; Clark et al. 1975). With the improvement of X-ray instruments, sources with 
X-ray luminosities lower than 10 34 ' 5 erg s _1 were discovered with the Einstein Observatory 
(Hertz and Grindlay 1983) in 1980s. The study of X-ray sources in Galactic globular clusters 
with the Chandra X-ray Observatory can reach luminosities down to Lx = 10 29 ~ 30 erg s _1 
in the 0.5-6 keV band. The most luminous X-ray sources in globular clusters are neutron 
stars accreting matter from a companion star, and some dim X-ray sources with luminosities 
Lx > 10 32 erg s _1 may be neutron stars in quiescence. According to previous research, dim 
X-ray sources with X-ray luminosities ranging from 10 30 erg s' 1 to 10 32 erg s" 1 are roughly 
identified as accreting white dwarfs or X-ray active low-mass main sequence star binaries 
(see e.g., Grindlay et al. 2001; Pooley et al. 2002, 2003; Bassa et al. 2004; Heinke et al. 
2005; Verbunt & Lewin 2006). 

The high stellar densities in globular cluster cores make secure optical identifications 
of X-ray sources difficult. However, the high spatial resolution of Chandra allows us to 
determine X-ray source positions to sub-arcsecond levels, which can be used to identify 
optical counterparts in Hubble Space Telescope (HST) observations of the densest regions 
of globular clusters. Detailed identifications of low X-ray luminosity sources by using multi- 
wavelength approach has been applied in some nearby globular clusters and some globular 
clusters with low central density (e.g. NGC 6397 by Grindlay et al. 2001, 47 Tuc by Heinke 
et al. 2005, M4 by Bassa et al. 2004, NGC 288 by Kong et al. 2006, NGC2808 by Servillat 
et al. 2008, M55 & NGC6366 by Bassa et al. 2008 and M12 by Lu et al. 2009). These 
studies suggest that the dim sources in globular clusters are quiescent low-mass X-ray binaries 
(qLMXBs; Wijnands et al. 2005 and references therein), cataclysmic variables (CVs; Cool 
et al. 1995; Bailyn et al. 1990), chromospherically active binaries (AB; Cool et al. 2002) or 
millisecond pulsars (MSPs; Grindlay et al. 2002). 

CVs consist of two components: a white dwarf and a main sequence or sub-giant star. 
Because of the accretion on the white dwarf, CVs are bluer than main sequence stars in the 
optical, and have higher Ha and ultraviolet emission. Quiescent LMXBs are binaries with a 
neutron star but the accretion rates are much lower than in persistent LMXB. The optical 
colors of qLMXBs are between CVs and main sequence. Chromospherically ABs have strong 
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coronal activity because of faster rotation periods. MSPs are neutron stars that emit radio 
emission with millisecond periods, and often have white dwarf or main sequence companions. 

The study of Verbunt(2002) suggested that CVs are formed from close encounters while 
ABs from primordial binaries. The study of Pooley et al.(2003; see also Heinke et al. 2003) 
presents conclusive evidence of the link between the number of close binaries observed in 
X-rays ( Lx(o.5-6fc e y) > 4 x 10 30 ergs s _1 ) and the encounter rate of the globular cluster. 
It also suggests that CVs are mostly formed from stellar encounters. Theoretical results of 
Davies (1997) also show that the progenitor binaries of CVs from wide primordial binaries 
are destroyed in the cores of globular clusters with high core densities (> 10 3 stars pc~ 3 ). 

Binary stars play important roles in star and globular cluster formation and evolution. 
The binding energy of these systems may be greater than that of the whole cluster. The 
dynamical process and the energy transfer of the binaries and the stars can strongly affect 
the dynamical evolution of the globular cluster (Hills 1975; Hills 1976; Pryor et al. 1991). 
Among clusters with a low core density and low encounter rate T, the number of faint X- 
ray sources found in NGC 288 (Kong et al. 2006) is higher than the prediction based on 
stellar encounter frequency (Verbunt 2003; Pooley et al. 2003), which has been speculated 
as the evidence for primordial binaries. Pooley and Hut (2006) describe the number of X-ray 
sources (N) in globular clusters using a combination of the collision number (r) and mass 
(M): N=ar+bM, where M is the cluster mass. The best solution for a predicted number of 
X-ray sources in a globular is given by a,— 1.2 and b=l.l (when T and M are normalized on 
the values for M4; Verbunt et al. 2007; Bassa et al. 2008). 

We report here on the Chandra and HST observations of the low core density globular 
clusters NGC6144 (p c = lO 2 - 23 L p C - 3 ) and E3 (p c = lO L11 L p C - 3 ; Harris, W.E. 1996); the 
Chandra observation and data analysis are described in §2 , the HST observation in §3 , 
sources identification in §4 and discussion in §5 . 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 

The galactic globular cluster NGC6144 is located at a J200 o=lQ h 27 m U.l s <5j 20 oo=-26 01'29' 
(Alcaino. 1980) and E3 is located at a J200 o=09 A 20 m 59.3 s 5/ 2 ooo=-77°16'57" (van den Bergh, 
Demers, and Kunkel 1980). Parameters of NGC6144 and E3 used in this paper are given in 
Table 1. 

NGC6144 was observed for 55410s on 2007 July 6, and E3 was observed for 20110s on 
2007 May 13, both with the Advanced CCD Imaging Spectrometer (ACIS) back-illuminated 
S3 chip on the Chandra X-ray Observatory. The field-of-view of the S3 chip (~ 8.3 x 8.3 
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arcmin 2 ) covers the entire half-mass radii of the two globular clusters. The data were collected 
with nominal frame time 3.2 sec. We only focus on the sources detected on the S3 chip. 

2.1. Data Reduction and Sources Detection 

The data analysis was done with Chandra Interactive Analysis of Observations (CIAO) 
3.4 software developed by the ChandraX-r&y Center. We used CIAO wavelet-based wavdetect 
(Freeman et al. 2002) tool to detect discrete sources in 0.3-7 keV band. The light curve of 
the whole ACIS-S1 chip indicates there are no background flares. We also performed source 
detection on the 0.3-1, 1-2, and 2-7 keV images for both NGC6144 and E3. The detection 
threshold was set to be 10~ 6 so that the false detection resulting from noise fluctuation is 
less than one from a 10 6 -pixel ACIS-S3 image. We set the wavelet radii from 1 to 16 pixels 
increasing by a factor of y/2. A total of 36 X-ray sources for NGC6144 and 16 X-ray sources 
for E3 were detected. Among the sources detected by Chandra , 6 lie within the half-mass 
radius(l'.62) of NGC6144 while 3 lie within the half-mass radius(2'.06) of E3 (Figure 1). The 
source counts were extracted from an elliptic region which covers 3a of the Chandra point 
spread function. The background count rate was extracted from the region without source 
outside the half-mass radii on the S3 chip. 

Table 2 lists the detected Chandra X-ray sources in NGC6144 and E3. Sources were 
named by the total counts in 0.3-7keV band, and 6 sources out of NGC6144 (CX01 to 
CX06) as well as 3 sources out of E3 (CX01, CX02, CX03) lie in their own half-mass radius, 
respectively. The columns give the sources' name, the coordinates (J2000.0), the net counts 
of three energy bands which are soft (0.3-lkeV), medium (l-2keV), hard (2-7keV), the 0.5-2.5 
and 0.3-7keV counts and unabsorbed fluxes. The vignetting caused by off-axis flux from all 
sources within the half-mass radius accounts for less than 3% of the change of flux. This 
value is less than the count error of each source. We therefore did not include vignetting 
calibration of flux of the source. We use a power-law spectrum with a photon index of 1.5 
and N H = 3.0 x 10 21 cm~ 2 for NGC6144 (estimating from Figure 2, see §2.2), photon index 
of 2 and N H = 9.8 x 10 20 cm" 2 for E3 (Kalberla et al. 2005) to estimate the unabsorbed 
fluxes. 

We estimate the completeness limit to be F x = 1-62 x 10~ 15 ergs s _1 cm~ 2 for NGC6144 
and F x = 2.65 x 10~ 15 ergs s~ 1 cm~ 2 for E3 in 0.3-7keV energy band. This corresponds to > 
10 counts and > 7 counts for NGC6144 and E3 anywhere on the S3 chip respectively. To 
obtain a rough estimate of the expected number of background sources expected to be found 
within the half-mass radius we proceed as follows. The area of the S3 chip that we analyze 
is 8.3' x 8.3'. Writing the half-mass radius in arcminutes as x, and the number of sources 
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detected outside the half-mass radius as N, we obtain the expected number of background 
sources within the half-mass radius as iV7rx 2 /(8.3 2 — rrx 2 ). This gives 4 and 3 background 
sources expected within the half-mass radii of NGC6144 and E3, respectively. The numbers 
of background sources estimated from the flux distribution in the Chandra Deep Field South 
(4 for NGC6144 and 2 for E3 respectively; Giacconi et al. 2001) are compatible with this. 
With 6 and 3 sources observed, was cannot exclude on the basis of the source numbers alone 
that all detected sources are due to background. 



2.2. Count Rates, Light Curves and Spectral fitting 

We only extracted the energy spectra and light curves for sources with more than 100 
counts: CX01, CX07, and CX08 in NGC6144. We rebin the spectra with 10 counts per 
spectral bin for CX01 and CX07, but 5 counts per spectral bin for CX08. Because of the 
low count numbers, we apply C-statistics in the spectral fitting (Nousek & Shue 1989). 
We fit these three sources with simple absorbed power-law and thermal bremsstrahlung 
model(Table 3). The fitting of each source with both models seems equally well from the 
C-statistics. The light curves of CX01 and CX07 are 2ks per bin and CX08 3ks per bin. We 
did not find significant variation for any of the sources from their light curves. 

The source counts are too few for the remaining sources to derive spectral parameters 
with meaningful results. Therefore, we also use hardness ratios to get rough X-ray spectral 
properties for all detected X-ray sources. These ratios were based on the source counts in 
three energy bands: S (0.3-1.0keV), M (l-2keV), and H (2-7keV). We analyzed two hardness 
ratios which are defined as HR1 = (M-S)/(M+S) and HR2=(H-S)/(H+S). Figure 2 shows the 
X-ray color-color diagram for NGC6144 (left) and E3 (right) and is overlaid with four lines 
showing the tracks of different Nh and photon index. From figure 2, a spectrum with photo 
index 1.5 and N H = 3.0 x 10 21 cm~ 2 would give more accurate unabsorbed flux than that 
with N H = 1.28 x 10 21 cm~ 2 given by Kalberla et al. (2005). The color-magnitude diagrams 
for all detected X-ray sources are displayed in Figure 3. 



3. OPTICAL OBSERVATIONS AND DATA ANALYSIS 

Optical observations for NGC6144 and E3 were taken with the HST Wide Field Camera 
of the Advanced Camera Surveys (ACS-WFC) with a field of view of 3.4' x 3.4', on 2006 
April 15 (Proposal ID 10775). Both globular clusters were observed in F606W (V 60 6) and 
F814 (J 814 ) filters. For NGC6144, the exposure time with the F606W and F814W filters are 
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1725s and 1775s. For E3, the exposure time with F606W and F814W filters are both 405s. 
Two X-ray sources in NGC6144 (CX03 and CX05) are also included in an archival HST 
Wide Field and Planetary Camera 2 (WFPC2) observation (Proposal ID 11014). F336W 
(C/ 3 3 6 ), F439W (£439), F675W (R 675 ), and F656N(Ha) filters were used in this observation. 
The exposure times were 3100s for F336W, 1700s for F439W, 3800s for F656N, and 900s for 
F675W filter respectively. Figure 4 shows the ACS-WFC field of view covering the entire 
half-mass radius of NGC6144 and most of the half-mass radius of E3, along with the field of 
view of the WFPC2 image for NGC6144 is also shown in Figure 4. 



3.1. Photometry 
ACS-WFC 

The photometry was done with the DOLPHOT package which is adapted from HSTphot 
(Dolphin 2000) for the use of ACS data. Because of the resampling of images, the drizzled 
images produce suboptimal photometry; therefore we used the individual flat-fielded images 
for photometric reference and combined images for astrometry. The acsmask masks out bad 
pixels from data quality images provided by the Space Telescope Science Institute (STScI), 
and we got images in units of electrons on the raw images. We employed splitgroups to 
split the ACS-WFC data with two chips into two extension files. The calcsky can create 
sky images for background determination. We use acsf itdistort to get the shift from the 
reference image to others. 

There were five and four exposures with different physical offsets per filter for NGC6144 
and E3, respectively. Two exposures with the same filter during the observations of NGC6144 
and E3 have no relative offset. We therefore use them as reference images. We gave the shift 
value for each chip, but the magnification and rotation are corrected by DOLPHOT. We 
ran photometry using DOLPHOT with ACS-WFC point-spread functions. DOLPHOT can 
analyse images within the same field and give individual photometry for each filter. Output 
photometry data were filtered to eliminate artifacts and obvious false detections. 

The photometries by Sarajedini et al. (2007) of NGC6144 are consistent with what 
we obtained by using DOLPHOT, but they do not provide magnitudes for all stars that are 
located in our error circles. However, in the photometry of E3 done by DOLPHOT, the stars 
were slightly distorted. The list from Sarajedini et al. provides better photometry for E3; 
we therefore use the result from Sarajedini et al. for the magnitudes of optical counterparts, 
whose color-magnitude diagrams (CMDs) are shown in Figure 5. 
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WFPC2 

Because the spatial resolution of ACS-WFC ( 0".05) is higher than WFPC2 ( 0".l for WF 
and 0".05 for PC), CX03a, CX03b, CX05a , and CX05b cannot be well resolved in WFPC2 
observations. We therefore provide the source list from DOLPHOT and use HSTPHOT to 
do WFPC2 photometry. The position of stars in the source list was calibrated by center in 
IRAF. Some steps were performed before applying HSTPHOT. We use mask for data quality 
image provided by the STScI to mask bad pixel. We then remove the cosmic ray on the 
images by using crmask. The sky value is determined by getsky. The command hotpixels 
can remove all hot pixels. The measurements were done using a predetermined star list. The 
CMDs of WFPC2 F336, F656N, and F675 for NGC6144 are presented in Figure 6. 

3.2. Astrometry 

In order to search for optical counterparts of Chandra X-ray sources in the fields, we 
have to improve the astrometry of each image. For NGC6144, we use relative astrometry 
instead of absolute astrometry to achieve a more accurate position. But for E3, the stellar 
density is relatively low. We used 2MASS (Two Micron All Sky Survey) point source catalog 
as reference to improve the astrometry. 

Astrometry for NGC61U 

For relative astrometry, we need to find sources observed both in the X-ray and optical 
bands, and determine the relative shift between the HST ACS-WFC image and Chandra 
image. The optical stellar density of NGC6144 in the ACS-WFC field of view is too high to 
find reference stars; we can identify CX28 with a reliable optical source only, which would not 
be accurate enough to determine precise astrometry. Instead, we use an image of NGC6144 
taken with a wide field imager. This image was taken by the European Southern Observatory 
(ESO) Max Planck Gesellschaft (MPI) 2.2m telescope with the Wide Field Imager (WFI) 
on 2002 February 21. The exposure time is 149.9s with the B-band filter. The images of 
the WFI were reduced using the Image Reduction and Analysis Facility (IRAF). The images 
were bias-subtracted and flat-fielded. To match coordinates on the wide field image with 
the HST image, we identify 23 isolated stars on the HST image to overlay with the same 
set stars on the WFI image. This process was done by ccmap in IRAF with a registration 
error of 0".06. There are three isolated optical sources in the WFI corresponding to Chandra 
sources CX21, CX31 and CX32. The difference of coordinates on the HST drizzled image 
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to the Chandra image corresponds to a shift of 0".156 in R.A. and a shift of 0".219 in decl 
with an error of 0".212. 

To search for possible optical counterparts, we combine three errors from astrometry and 
source detections to determine the radius of the error circle: the registration error between 
the HST image and the WFI image, the error resulted from astrometry difference between 
Chandra and WFI, and the uncertainty of X-ray source positions detected by CIAO. The final 
uncertainty on the position of X-ray sources is the quadratic sum of the above-mentioned 
uncertainties of the X-ray sources. The final astrometric solution gives an error of 0".22. 
The following identification will be focused on sources within the 95% error circles. 

Astrometry for E3 

To improve the absolute pointing accuracy of X-ray and optical frame, we use 2MASS 
All-Sky Catalog for Point Sources as reference for absolute position on both optical and X-ray 
images. We found 66 point sources on the HST ACS-WFC image in a 3 x 3 arcmin 2 region 
around the center of E3 that matched entries in the 2MASS catalog. This gives positional 
rms residuals of 0".173 in R.A. and 0".155 in decl. For the ACS-WFC image, the astrometric 
solution using the 66 stars yields rms residuals of 0".134 in R.A. and 0".162 in decl. 

For the Chandra image, we found two X-ray sources (E3-CX01 and E3-CX14) match the 
2MASS catalog. The calibration (using E3-CX01, E3-CX14 and the 2MASS counterparts) 
of coordinates on the Chandra image corresponds to a shift of 0".134 in R.A. and a shift of 
0".575 in decl with residuals of 0".062 in R.A. and 0".425 in decl. 

Combining the two 2MASS stars with rms residuals of 0".183 in R.A. and 0".187 in 
decl. gives the final astrometric solution with residuals of 0".487 in R.A. and 0".560 in decl. 
For the astrometric calibration of X-ray and optical images, we get a 0".805 positional error. 
The error is the quadratic sum of the uncertainty of 2MASS sources for astrometry and 
the absolute astrometry for HST and Chandra and the uncertainty of X-ray source position 
detected by CIAO. 

We later identify E3-CX01 in Chandra image with an HST object with certainty. This 
source in both images shares one same 2MASS star, so the uncertainty of the position of 
E3-CX01 is simply that given by wavdetect. However, while one optical source has been 
identified as a confirmed HST counterpart to a Chandra source (E3-CX01), the error circles 
for other Chandra sources (E3-CX02 and E3 CX03) on the HST image can be reduced to the 
sum of the wavdetect error for that source and the wavdetect error for the E3-CX01 (since 
relative HST positional errors are negligible, typically < 0".01). Using the wavdetect errors 
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in Table 2, we find the positional error for E3-CX02 is 0".43 and E3-CX03 0".38. These 
errors are less than those estimated from absolute astrometry. Figure 7 gives sources within 
the 95% error circles. 



4. SOURCES IDENTIFICATION 

To identify and classify sources in NGC6144 and E3, we analyse data from both X-ray 
and optical observations. We search for optical counterparts within the 95% Chandra error 
circle. We discuss all X-ray sources within the ACS-WFC field of view but only consider 
sources within the half-mass radius to be cluster members in §5 . For multiple HST sources 
within the Chandra error circles, all the candidates in 95% error circles are summarized 
in Table 4, Table 5 provide WFPC2 photometry for the candidates in the error circles of 
CX03 and CX05, and the finding chart is shown in Figure 7. Star just above the main 
sequence could be binaries in the cluster; stars left of the main-sequence could be binaries 
with accretion disks or hot white dwarfs. Other possibilities are foreground or background 
stars and background galaxies or AGNs. 

The spectral characteristic is described by the hardness ratios (soft, medium and hard) 
in the X-ray color-color diagrams and color-magnitude diagrams. We also made a plot of the 
X-ray luminosity with respect to the optical V-band absolute magnitude (Figure 8). Based 
on previous studies on M4 (Bassa et al. 2004), NGC288 (Kong et al. 2006) and M55 & 
NGC6366 (Bassa et al. 2008). The left diagonal line in Figure 8 roughly separates CVs and 
ABs. We calculated the X-ray flux to optical V-band apparent magnitudes as: log (fx/ fv) 
= log fx + 5.67 +0.4w 6 o6 (Green et al. 2004); fx is derived in the 0.5-2.5 keV band; we can 
roughly classify some sources with fx/ fv > 0.8 (Cool et al. 1995) and My < 11 (Patterson 
1984) as CVs. 

The Chandra X-ray sources NGC6144-CX28 and E3-CX03 are extended objects in the 
optical images and the optical color-magnitude diagrams also show that they are redder 
than main sequence stars, indicating they are background galaxies. The X-ray luminosity 
of NGC6144-CX01 (£x(o.5-6fc e y) ~ 1-8 x 10 32 erg s -1 ) is too high for a binary (typically 
L X (o.5-6keV) < 10 30 erg s _1 from Verbunt & Lewin 2006) of two main sequence stars; in ad- 
dition, it has a high X-ray to optical luminosity ratio and a constant light curve. The optical 
counterpart to CX01 might be a chance coincidence with a background AGN. However, even 
in the case when we decrease the optical magnitude so the counterpart becomes fainter and 
CX01 shifts to the left in Figure 8, CX01 is still far from the parameter space occupied by 
galaxies in previous study. We suggest that NGC6144-CX01 is a possible CV, even though 
the optical counterpart on CMD lies on the main sequence. The spectrum of NGC6144- 



- 10 - 



CX02 is relatively hard in X-rays. CX02 is bluer in optical and has a high X-ray to optical 
flux ratio. The optical magnitude suggests it is possibly a faint CV or a galaxy judging from 
its position in Figure 5 (similar to position in the CMD with M55-CX2 from Bassa et al. 
2008). There are two optical counterparts in the error circle of NGC6144-CX03, CX03a and 
CX03b. The one possible optical counterpart (CX03a) is redder than the main sequence and 
the ratios of X-ray flux to optical V-band magnitudes is low (<0.2). However, CX03a is un- 
likely to be a point source, CX03a is a saturated star and may contain more stars. Certainly 
identification would not be possible without further observations. The other possible optical 
counterpart (CX03b) is bluer than the main sequence with Ha emission, even the ratios of 
X-ray flux to optical V-band magnitudes is low (<0.5). As a result, we suggest that CX03 
might be classified as a CV. 

The X-ray to optical flux ratio of NGC6144-CX04 is lower than CX03, and the optical 
counterparts (CX04a and CX04b) also lie on the main sequence; therefore it is also an 
possible AB. There are three possible optical counterparts in the error circle of the faint 
source NGC6144-CX05. The brightest counterpart (CX05c) saturated the HST ACS-WFC 
detector; we can only analyse the other two counterparts. Two stars in the error circle 
of CX05, CX05a and CX05b, are on the main sequence. The X-ray to optical flux ratio 
for CX05a is 0.66; it may be an active binary. The X-ray to optical flux ratio of CX05b is 
somewhat high for an active binary. In the photometry of WFPC2 image, CX05a and CX05b 
are too faint (S/N < 2) to provide reliable magnitudes. Classification of CX05b wouldn't 
be possible without further information. The lower limit (magnitude) of the saturated star 
(CX05c) in the error circle was 17.96 for C/ 336 , 17.51 for £439, 16.24 for Ha, 17.76 for R 656 , 
18.84 for V606, and 15.63 for 7gi4- The I-band is three magnitudes brighter than V-band; 
hence it looks like a red giant. However, since it is saturated on the detector, we cannot 
assign it a correct color. 

For NGC6144-CX07, the most luminous X-ray source in NGC6144 has two possible 
optical counterparts: CX07a and CX07b. The optical color-magnitude diagram of CX07a 
is bluer and CX07b redder, and the X-ray flux to optical flux ratio of CX07a and CX07b is 
3.4 and 608.9 respectively. We identify CX07b as an AGN. Based on the optical results and 
X-ray hardness ratios of NGC6144-CX07, it could be either a CV or an AGN with CX07a 
as the optical counterpart. Certainly identification would not be possible without further 
observations. We found two possible optical counterparts for NGC6144-CX08: CX08a and 
CX08b. Both of them appear on the main sequence, but the X-ray flux is still too high for 
CX08 to be an AB. The X-ray to optical flux ratio also indicates that NGC6144-CX08 is 
more likely to be a CV. 

Based on the optical source density near the X-ray error circles the possibility of at 
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least one the chance of at least one star falling by chance into the error circle of previous 
identification of ABs in this Section are: 35% for CX03, 45% for CX04, 91% for CX05, and 
44% for CX08. We therefore cannot exclude the possibility that all our ABs are positional 
coincidence. Indeed, the upper limit for ABs (the dashed line in Figure 8) given by Bassa et 
al. (2008; Equation 3) indicates that CX04a is more likely to be a possible AB than other 
sources higher than this limit. This limit for ABs is based on previous studies for stars in 
the solar neighborhood. We take this limit for AB in NGC6144 because the density of the 
cluster is low and the identified ABs in low core density globular cluster M12 (CX4b; Lu et 
al. 2009), M55 (CX7), and NGC6366 (CX4 and CX5; Bassa et al. 2008) are also lower than 
this limit. However, note that previous studies of ABs in high core density globular clusters 
can be higher than this limit. The position of CX04a in the CMD appears to be a subgiant, 
similar to the position of M55-CX7 in the CMD from Bassa et al. (2008). Some X-ray error 
circles contain more than one optical object without unusual features. They clearly require 
at least one optical object must be a coincidence. 

The X-ray color-color diagram shows that E3-CX02 has a relatively hard X-ray spectrum 
compared to E3-CX01 and E3-CX03. The optical image of E3-CX02 also shows that it is a 
blue star and the high X-ray flux to optical magnitude suggests this is a possible CV. But 
the color of E3-02 is bluer than turnoff. And it is one magnitude bluer than main sequence 
at same I-band magnitude. It appears bluer compared to CVs previously studied in globular 
clusters. Certainly identification would not be possible without further observations. 

There are two unclassified X-ray sources (CX06 and CX30) in the ACS-WFC field of 
view of NGC6144 that no optical counterparts are found. The optical counterpart for E3- 
CX01 is too bright and saturates the CCD, it may be a foreground star. All other sources 
can not be classified without further observation. These two unidentified sources (CX06 and 
CX30) may be background sources which are too faint to be detected by HST ACS-WFC in 
the W606F and W814F filters. 

We shift the error circle on the optical image of NGC6144 and E3 to eight different 
directions (respectively) about 5" for NGC6144 and 10" for E3 so that every time some 
sources would lie in the circle. According to the average number of sources found in shifting 
the circle (average over 8 trials in 8 directions), Table 6 listed the probabilities that the 
sources in each error circle come from positional coincidence. For CX04, based on the low 
probability (<6 %) of finding two or more sources within the error circle, one of the two 
sources (CX04a, CX04b) is probably the real counterpart. For CX03 (the identified CV in 
this section), the probability that all the optical sources within the error circle are chance 
coincidence is also relatively low (<10%). CX02 and CX05 are two marginal cases and we 
cannot exclude that all sources within the error circle are not real counterpart. Moreover, 
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the probabilities that all sources in our error circle are coincidence (not true counterpart) are 
0.06% for NGC6144 and 0.92% for E3. Those low probabilities suggest that the potential 
candidates of optical counterparts in our error circle are very likely to be true counterparts. 

5. Discussion 

Based on the identification in §4 and the probability of positional coincidence in Table 
6, there is one good CV candidate (CX03) and two possible CVs (CX01 and CX02) within 
the half-mass radius of NGC6144. We also found one source (CX04) as active binary. The 
identification in §4 indicates that at least 2 sources (CX03, CX04) are probable members 
associated with NGC6144 while other 2 sources (CX02, CX05) are marginal. We therefore 
assume 3 sources as cluster members of NGC6144 in the following discussion. If we take 
3.5 sources to be chance coincidence in 6 error circles (from CX01 to CX06, see Table 6), 
the probability that at least 3 sources are cluster members is higher than 86%. Three 
sources in globular cluster E3 were detected within the half-mass radius, with one insecure 
CV. The study by Pooley et al. (2003) indicates the relation between the stellar encounter 
rate and the number of X-ray sources with L X (o.5-6keV) > 4 x 10 30 ergs s -1 in globular 
clusters. The encounter rates are proportional to p c 2r c v rm,s where p c is the central density 
of globular cluster, r c the core radius, and v rms the central velocity dispersion. We estimate 
the encounter rate by Toe p l c ^r 2 c via v rms oc (p c r c)° 5 (Verbunt 2003). The encounter rate for 
NGC6144 (r 6 i44 = 0.07) is about 1/330 of 47 Tuc and 1/15 of M4. For E3 (T E3 = 0.001) it is 
about 1/15000 of 47 Tuc and 1/700 of M4. The relevant detected numbers of X-ray sources 
( L X (o.b-6keV) > 4 x 10 30 ergs s" 1 ) in 47 Tuc and M4 are 45 and 6 respectively (Pooley et 
al. 2003). It is clear that the X-ray sources number in low core density globular cluster, 
NGC6144 (2-3 sources) and E3 (one doubtful source), does not scale with T. 

We consider the number of X-ray sources with the collision number and mass (Pooley 
& Hut 2006). The best fit for a predicted number of X-ray sources in a globular cluster by 
Bassa et al. (2008) is N = 1.2r + l.lM h . Values for T and M h are normalized to the value 
of M4. According to this fitting, one X-ray source is predicted in NGC6144. 

We also repeated the fitting based on Bassa et al.(2008) with the same clusters (Pooley 
et al. 2003) as well as NGC288(^=8), M55(^=7), NGC6366( / u fe =4), NGC6144(/if,=4.5), 
and NGC6218(/i fe =4; Lu et al. 2009) where p b is the expected number of background sources 
(Verbunt et al. 2007). Verbunt et al.(2007) suggest a Poisson distribution to overcome the 
small-number statistics of clusters with small source numbers. According to the Poisson 
distribution, the probability of observing N X-ray sources when p is expected is P(N, p)=p N 
e~ M /N!. The fitting procedure consists of varying a and b in N = aF + bM^ to maximize 



- 13 - 



the probability P=Y[j[P(N c ,fj, c )P(N b ,fj, b )]j ,where j indexes the clusters, N b the number of 
background sources not related to the cluster, N c the number of cluster sources within the 
half- mass radius (subtract N b from observed number of X-ray sources), and ji c the expected 
number of cluster sources. For these fittings we assumed the lower limit of source number 
to be 2 for NGC288, 3 for M55, 2 for NGC6366, 2 for NGC6218 (Lu et al. 2009) and 3 for 
NGC6144. The lower limit of the other clusters sources is 0. The best solution is found for 
N = 1.23r + 1.14M ft with the total probability P=e~ 57 (see Verbunt et al. 2007). This result 
is consistent with Bassa et al. (2008) even when we include NGC6144 and NGC6218. Figure 
9 gives the comparison between the observed number of sources (background-subtracted) 
and the best fit number. 

We fitted the data anew, excluding NGC6218, NGC288, M55 and NGC6366. The best 
solution is found to be iV = 1.41T + O.Q2M h with P=e~ 37 (without lower limit to the source 
number of NGC6144) or iV = 1.32r + 0.83M h with P=e~ 38 (with lower limit to the source 
number of NGC6144 equaling to 3). According to the previous fits result N = 1.23T+l.lAM h , 
the first term, which is the stellar encounter rate, yields 0.1 expected sources (r6i44 = 0.07), 
the second term, which is the mass of NGC6144, yields 1.0 expected source (M6144 = 0.87). 
The dependence of the X-ray sources in NGC6144 on mass is significantly shown. 

If E3-CX02 is the member of E3, then the number of identified X-ray sources in globular 
cluster E3 is still higher than the prediction based on the relation between the mass of the 
cluster and the number of sources. This may be due to the structure of the cluster is now 
being disrupted under the galactic tidal force (van den Bergh et al. 1980). The encounter 
rate and the mass of E3 are 0.001 and 0.022 (normalized with M4), respectively. The X-ray 
sources in E3 should have more dependence on mass than T. 

The scaling of the number of detected source number found with £v(o.5-6fceV) > 4 x 10 30 
ergs s _1 (Pooley et al. 2003; see also Verbunt et al. 2002; Heinke et al. 2003) suggests 
that CVs are mostly made from close encounters. The population of CVs may be caused by 
either tidal capture or interactions involving primordial binaries (Davies 1997). In a low core 
density globular cluster, primordial CV-progenitor binaries are likely to survive long enough 
to form CVs in the cluster core. NGC6144 and E3 are the two clusters selected to study the 
faint X-ray sources population in low core density globular clusters. By combining both the 
Chandra and i^STdata, we found possible CVs in NGC 6144, and no good CV candidates in 
E3. Since these CVs are unlikely to be formed via dynamical processes (stellar encounter), 
they must be evolved from primordial binaries as predicted by Davies (1997). 

There is evidence that the number of X-ray sources in low core density globular cluster 
does not scale with V are also shown in the globular cluster NGC288 (also Kong et al. 2006). 
The study of M55 and NGC6366 by Bassa et al. (2008) clarifies the dependence of the 
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number of X-ray sources not only on T but also on the mass of the globular cluster. In brief, 
there are evidences that close binaries in globular cluster could be primordial in origin. 
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grateful to Andrew Dolphin for his DOLPHOT code and some assistance to solve photomet- 
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Inc., under NASA contract NAS5-26555. These observations are associated with HST pro- 
posal number 10775 and 11014. This project is supported by the National Science Council 
of Republic of China, Taiwan through grant NSC 96-21 12-M-007-037-MY3. 
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Fig. I.— Chandra X-ray sources on the entire ACIS-S3 chips with 0.3-7 keV of NGC6144 
(left) and E3 (right). The square is the field of view of HST ACS-WFC. The large circle is 
the half-mass radius of the cluster and the smaller circle the core. All detected X-ray sources 
are marked with source number. 
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Fig. 2. — X-ray color-color diagram for most of the X-ray sources detected in Chandra S3 
chip of NGC6144 (left) and E3 (right). The sources marked with error bars and source 
names are the ones that lie in the field of view of HST ACS-WFC. We fit the data with 
power law model. The values of N H in NGC6144 are 1.5 x 10 21 , 3.0 x 10 21 , 5.0 x 10 21 , and 
6.0 x 10 21 cm~ 2 (from left to right). We also applied power-low model on E3; the values of 
N H are 1.5 x 10 21 , 3.0 x 10 21 , and 6.0 x 10 21 (from left to right). For each power-law model, 
the photon index varies as (from bottom to top) a=3.0, 2.75, 2.5, 2.25, 2.0, 1.75, 1.5, 1.25. 
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Fig. 3. — X-ray color-magnitude diagram for most of the X-ray sources detected in Chandra 
S3 chip of NGC6144 (left) and E3 (right). The X-ray color is defined as the logarithm of the 
ratio of 0.3-1.0 keV (X so / t ) counts to 2.0-7 keV (X hard ) counts, and the magnitude is the 
logarithm of 1.0-2.0 keV (X mid ) counts. The sources marked with error bars are the ones 
that lie in the field of view of HST ACS-WFC. 



NGC6144 E3 




Fig. 4. — The HST ACS-WFC images with X-ray sources detected by Chandra with 0.3-7 
keV of NGC6144 (left) and E3 (right). The square is the field of view of HST ACS-WFC. 
The polygon on NGC6144 is the field of view of HST WFPC2. Half-mass and core radii are 
shown. 
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Fig. 5. — Optical color-magnitude of the HST ACS-WFC observation of NGC6144(left) and 
E3 (right). Possible X-ray source candidates for optical counterparts are labeled with source 
name. The letters after source names refer to the candidate optical counterparts to the X-ray 
sources. 
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Fig. 6.— Optical color-magnitude of the HST ACS-WFPC2 observation of NGC6144. Only 
CX03a and CX03b provide reliable magnitudes. 
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Fig. 7. — The (5x5 arcsec 2 ) view the candidate optical counterparts of NGC6144 and 
(10 x 10 arcsec 2 ) view of E3. These images were taken from HST ACS-WFC observations in 
F606W. The 95% error circles were overlaid for the Chandra source position. All images of 
NGC6144 use the same gray scale throughout except for CX03. The images of E3 use the 
same gray scale throughout but different from that used for NGC6144. 
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Fig. 8. — X-ray luminosity to V-band absolute magnitude diagram. The X-ray luminosity 
L x is in the range from 0.5keV to 2.5keV. The classified X-ray sources from NGC288 , 47Tuc, 
NGC6121, and NGC6752 are also plotted with NGC6144 and E3. Optical magnitude of V 
band for NGC6144 and E3 are observed using HST F606W filter. The solid line was used in 
earlier papers (e.g. Bassa et al. 2004) to separate cataclysmic variables above it from active 
binaries below it. However, the maximum X-ray luminosity observed for active binaries near 
the Sun is about 0.001 L bol (Verbunt et al . 2007, Bassa et al. 2008), and this line is indicated 
with a dashed line (based on L bo i from the 11.2 Gyr isochrone for metal-poor z=0.001 stars 
as computed by Girardi et al. 2000). 
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Fig. 9. — Comparison between the observed number of sources (symbols: □) and the best 
fit number (symbols: x). The position of each cluster is given by T and Mass. The scales 
for the number of sources are on the top left of the figure. The best fit gives number of each 
cluster by iV = 1.2T + l.lM h . The size of the symbols is proportional to their number. The 
mass and encounter number(r) of clusters are normalized with M4(NGC6121). 
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Table 1: Parameters of NGC6144 and E3 



Cluster 


r c 


r h 


d (kpc) 


logp c 


Eb-v 


N H cm.- 2 * 


NGC6144 


0.94' 


1.62' 


8.5 


2.23 


0.36 


1.28 x 10 21 


E3 


1.87' 


2.06' 


4.3 


1.11 


0.3 


9.36 x 10 20 



NOTES. — Values of angular core radius r c , angular half-mass radius r h , 
distance d, luminosity densities p c , reddening Eb-v, Neutral hydrogen 
column density Nh ■ 

* The interstellar column density of neutral hydrogen was taken by the 
average 21-cm emission from 1 degree 2 around the center of the globular 
cluster taken from The Leiden/ Argentine/Bonn (LAB) Survey of Galactic 
HI. 
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Table 2: Chandra Source Properties 
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NOTES. — List of Chandra X-ray sources in our observations of NGC6144 (upper part) and E3(bottom part). The 
positions of stars were given by event-2 file from Chandra observations and calibrated by using Two Micron All Sky 
Survey (2MASS) of Point Sources as reference. The position errors given by wavdetect are listed between brackets, 
in arcsec. The signal to noise ratios of photon counts were around 2-21. The source counts have been corrected for 
background. The first 6 sources are located within the half-mass radius of NGC6144 and are numbered according to total 
counts in 0.3-7keV band. The first 3 sources are located within the half-mass radius of E3 and are numbered according 
to total counts in 0.3-7keV band. The remaining sources are located outside the half-mass radius and are numbered 
according to total count in 0.3-7keV band. The unabsorbed flux is in units of 10 -15 ergs cm~ 2 s -1 . 
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Table 3: Spectral Fits of the Brightest Sources 
Source Name Model a kT/a Goodness / .3-7 d 



CX01 


PL 


0.781^ 


1.23-1.63 


73% 


5.2 




TB 


0.38lS;g 


< 49.84 


79% 


4.5 


CX07 


PL 


2.3418;™ 


1.70-2.20 


15% 


8.3 




TB 


i ci +0.51 
i - d± -0.39 


4.04-7.74 


10% 


4.4 


CX08 


PL 


4 16+ 2 - 27 

^• iu -0.34 


1.79-2.35 


78% 


3.7 




TB 


2 4n +0 - 08 

z '^ u -0.12 


2.93-7.89 


60% 


2.04 



NOTES. — The uncertainties of N H and kT/a are the 90% confidence. 
a PL: power-law.; TB: thermal bremsstrahlung. 
b in units of 10 21 cm -2 

c The goodness represents the fraction of simulations with a lower C-statistic 
than the data. If the spectrum and the simulations were produced by the 
same model, there is an equal probability of obtaining any given percentage. 
If this fraction is very large (perhaps 95%, or 99%) one may conclude that 
it is unlikely that the data is drawn from this model. 
d 0.3-7 keV unabsorbed flux in units of 10" 14 ergscm~ 2 s -1 . 
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Table 4: Optical Counterparts to Chandra X-ray Sources 



Source Name 


^606 


-^814 


fx/fv a 


Classification 13 


CX01 


23.25±0.01 


21.59±0.01 


27.81 


CV? 


CX02 


25.66±0.03 


24.38±0.03 


23.61 


CV? AGN? 


CX03a 


18.84±0.02 


17.39±0.01 


0.01 


AB? 


CX03b 


21.95±0.01 


20.68±0.01 


0.49 


CV 


CX04a 


18.69±0.01 


17.25±0.01 


0.01 


AB 


CX04b 


22.93±0.01 


21.33±0.01 


0.73 


? 


CX05a 


24.14±0.01 


22.30±0.01 


0.66 


? 


CX05b 


25.56±0.03 


23.41±0.02 


2.41 


? 


CX05c 


< 18.84 


< 15.63 


< 0.01 


AB? 


CX07a 


20.53±0.01 


19.66±0.01 


3.40 


CV7AGN? 


CX07b 


26.16±0.04 


23.11±0.01 


608.85 


AGN 


CX08a 


21.35±0.05 


20.05±0.08 


2.457 


? 


CX08b 


23.28±0.12 


21.44±0.12 


14.40 


CV? 


CX28 


24.43±0.016 


22.00±0.009 


2.11 


AGN 


E3-CX02 


22.79±0.03 


22.21±0.02 


2.37 


CV? AGN? 


E3-CX03 


25.06±0.09 


23.09±0.03 


17.02 


AGN 



a Ratio of X-ray to optical (Visoe) nux J using log(fxffv) = log fx + 5.67 + 0.4V606 (Green ct al. 2004); 
fx is derived in the 0.3-7 keV band. 

b CV: cataclysmic variable; AB: chromospherically active binary; AGN: active galactic nuclei. We consider 
all optical sources inside X-ray error circle as counterparts. 



Table 5: Optical Counterparts Observed by HST WFPC2 



Source Name 


^336 


-B439 


Ha 


-^675 


CX03a 


19.52±0.02 


19.23±0.02 


17.80±0.02 


18.35±0.02 


CX03b 


22.81±0.12 


22.74±0.13 


21.13±0.11 


21.45±0.10 


CX05a 


> 24 


> 24 


> 23 


> 23 


CX05b 


> 25 


> 25 


> 24 


> 24 


CX05c 


< 17.96 


< 17.51 


< 16.24 


< 17.76 



NOTES. - 



CX05a and CX05b are too faint to provide reliable magnitude. 
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Table 6: The Probability of Positional Coincidence 



Source Name 


Average Number a 


Probability b (%) 


CX01 


0.375 


31.27 


CX02 


0.375 


31.27 


CX03 


0.5 


9.00 


CX04 


0.375 


5.67 


CX05 


1.25 


15.49 


CX06 


0.625 




CX07 


0.25 


2.65 


CX08 


0.125 


0.89 


CX28 


0.5 


43.12 


CX30 


0.75 




E3-CX01 


0.125 


11.76 


E3-CX02 


0.125 


11.76 


E3-CX03 


0.125 


11.76 



a The error circle on the HST images of NGC6144 and E3 are shifted to eight directions (5" for NGC6144 
and 10" for E3). Here we list the average number of sources found in the circle over 8 trials. 
b The probability of finding the observed number (or more) sources inside the error circle, assuming 
Poisson distribution with a mean given by the average number. 



